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The problem of non-linear transport near a quantum phase transition is solved within the Landau 
theory for the dissipative insulator-superconductor phase transition in two dimensions. Using the non- 
equilibrium Schwinger round-trip Green function formalism, we obtain the scaling function for the non- 
linear conductivity in the quantum disordered regime. We find that the conductivity scales as E 2 at low 
field but crosses over at large fields to a universal constant on the order of e 2 /h. The crossover between 
these two regimes obtains when the length scale for the quantum fluctuations becomes comparable to 
that of the electric field within logarithmic accuracy. 



When an electric field is applied to a system near 
a quantum critical point, non-linear transport obtains. 
The non-linearity arises because the length scale [1] as- 
sociated with the electric field scales as Ie °c E^ 1 ^ 1+z \ 
where z is the dynamical exponent. As it is the ratio of 
the electric field length scale to the correlation length, 
£ w A - ", that enters the resultant DC resistivity, 

p(T, E)=f (A/T 1 /^, A/E 1 /^ 1 ^ , (1) 

a non-linear electrical response is inescapable. Here A 
is the distance from the critical point. Non-linear trans- 
port according to Eq. (1) is expected to hold as long as 
the temperature is low enough so that the length scale 
associated [7] with temperature, It ~ l/T 1 / 2 , exceeds 
that of the electric field, l/T 1 / 2 » Ie- Because of the 
additional factor of z that enters the electric field scaling 
contribution to the conductivity, simultaneous scaling of 
the resistivity with respect to temperature and electric 
field enables a direct determination of both the correla- 
tion and dynamical exponents, v and z, respectively and 
hence a complete characterization of the critical proper- 
ties. Despite its obvious importance, electric field scaling 
persists as the outstanding problem in quantum critical- 
ity because no theoretical account has been put forth 
to explain how quantum fluctuations conspire to yield 
non-linear transport. The primary theoretical hurdle is 
simple: A successful theory of non-linear transport must 
lie outside the standard Kubo/linear response formalism. 
Experimentally, the problem is complicated by the fact 
that in a wide range of systems exhibiting quantum criti- 
cal points [3-6], the I—V characteristics which are highly 
non-linear at small fields all become linear at large fields. 
In fact, the I — V characteristics for varying values of |A| 
all attain a universal slope as E — > oo. The vanishing 
of the non-linearity at high values of E implies that the 
scaling function for the resistivity is a highly non-trivial 
and non-monotonic function of the electric field. Hence, 
a successful theory of non-linear transport must uncloak 
how such non-monotonicity arises. 

Using a non-equilibrium formalism, we calculate ex- 



plicitly the scaling function for the non-linear conduc- 
tivity for both the large and small electric field limits. 
While the formalism can be applied to any critical the- 
ory, we focus on the Landau theory for the insulator- 
superconductor transition in which the dynamics are de- 
termined by an Ohmic bath; that is, z = 2. Within the 
Schwinger [8] round-trip double-time Green function for- 
malism (also known as the Kadanoff-Baym [9], Keldysh 
formalism [10]), we evaluate the current for arbitrary 
electric field simply by calculating the appropriate Green 
function. When temperature is the smallest parameter 
(quantum disordered regime) , a crucial ratio which enters 
both the conductivity and the inverse correlation length, 
to 2 , is the ratio 

Q = (e*£0 2/3 /A (2) 

(e* = 2e- the charge of a Cooper pair). For 
Q <C 1, corrections due to field are sub-dominant 
and to 2 = A/ln(l/A)) whereas for large field, m 2 = 
T(e* J B) 2 / 3 /ln(l/(e* J B) 2 / 3 , T a constant. In the quan- 
tum disordered regime, the static conductivity scales as 
Q 2 for Q <C 1, whereas in the large-field limit, Q 
1, linear transport emerges with a universal constant, 
u(E) = 0.46e 2 //i, that is in qualitative agreement with 
experiments [3-6]. 

The starting point for our analysis is the minimal 
Ginsburg-Landau action (in the imaginary time) 

F[4>] = J d 2 r J dri\ (v A(r, t)) ^(r,r)| 2 + 

|a^(r, r)| 2 + S Mr, r)| 2 + ((7/2) |^(r, r)| 4 } + i dis (3) 

required to model quantum fluctuations and dissipa- 
tion near the zero-resistance quantum critical point. In 
Eq. (3), A(r,r) is the vector potential, e* = 2e, 5 
is the bare distance to the quantum critical point. In 
terms of the Matsubara frequencies, the dissipation term, 
^dis = vJ2k,u n l w n||^(k, u n )\ 2 , corresponds to the phe- 
nomenological Ohmic model introduced by Caldeira and 
Leggett [11], in which r\ measures the strength of the dis- 
sipation. The order parameter ip(r,r) is the standard 
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two-component complex field whose expectation value is 
non-zero in the superconducting phase. Our motivation 
for introducing dissipation, that is expected to be strong 
(77 ~ 1), is two-fold. First, the second-derivative term 
now becomes irrelevant. Consequently, the zero temper- 
ature transition falls into the z — 2 universality class. 
Such a universality class describes a 2D insulator to high- 
T c superconductor transition in which superconductiv- 
ity is destroyed by impurities [12]. Second, as we will 
see, z = 2 dynamics are inherently easier to solve using 
the Schwinger [10,13] technique than is the correspond- 
ing z = 1 problem. Nonetheless, the generality of our 
conclusions leads us to believe that similar results must 
hold for the z = 1 case as well. 

As the problem we wish to treat is inherently out 
of equilibrium, we resort to the real-time 'round-trip' 
Schwinger formalism [8-10,13] which is ideally suited for 
solving problems in which an asymmetry exists between 
forward and backward evolution in time. We will be 
brief in our presentation of this technique as it is well- 
documented in the literature [10,13,14] and our notation 
will follow that of Lifshitz and Pitaevskii [13]. To use the 
real-time formalism, one needs to rewrite the action, Eq. 
(3, as an integral over the Keldysh contour. A constant 
in space and time electric field is assumed to enter the 
action via the vector potential, A(t) = — Ei. This gauge 
choice allows for all functions to be Fourier expandable 
in space, since full translational invariance is retained. 
We define the time ordered and anti-time ordered Green 
functions 

iG--(h,t 2 ) = FTl?\${?,t x \W M\) 
iG;+(t u t 2 ) = FT{[f^(vM),r(v'M)]), (4) 

where T and T are the time ordering and inverse time or- 
dering operators. Here FT represents the Fourier trans- 
form with respect to space. We will need also the non- 
time ordered functions 

iG-+(t 1 ,t 2 ) = FT(r(r',t 2 )iP(r,t 1 )) 
iG+-(h,t 2 ) = FT(iP(r,h)r(r',t 2 )), (5) 

the last of which is directly related to the expectation 
value of the current operator through 

The average values in these expressions are over all states 
of the system, not simply equilibrium ones. Since the four 
functions in Eqs. (4) and (5) forming a 2 by 2 matrix G, 
are not independent, one employs the so-called Keldysh 
rotation to work with the three independent functions 
G R (G A ) = G— - G-+(G+-),G K = G~+ + G+". The 
matrix of the exact Green functions G is connected to 
that of the Green functions for free quasi- particles G^ 



via the Dyson equation [13], G = G<°) +J J G^EG, in 
which the integration over internal time and space argu- 
ments is assumed and the self-energy matrix E is itself 
in general a complicated functional of G. For the static 
field, the resulting Green functions G + ~, G h , G K , taken 
at t = t\ = t 2 do not depend on t, reflecting the fact that 
the system is time translationally invariant as well. 

The non-linear conductivity is defined as the constant 
of proportionality between the current and the electric 
field. We will orient the field along x-axis and hence 
Jxx(E) = a xx (E)E. In the Schwinger formalism, the 
direct relationship [14] between G~ + and the retarded 
and advanced Green functions 

G p + (h,t 2 ) = - f dt s dt i G*(t u t 3 )X-+(t 3 ,t 4 )G£(t i ,t 2 ) (7) 

is obtained by formally solving the Dyson equation. In 
Eq. (7) , note the minus sign that is absent in the fermion 
problem [14]. Analogous expressions hold for G H (G K ) 
with E~+ replaced with E" 1 (E K ). In general, the self- 
energy E arises from the interaction U, £y, and the cou- 
pling to the dissipative bath E^. However, treating the 
interactions in the largc-N (mean- field) limit, we obtain 
that this approximation gives rise only to the renormal- 
ization of the bare distance to the critical point in the 
action [15], 

so that a new Gaussian action with the excitation spec- 
trum e(p) = p 2 + m 2 obtains. There is no con- 
tribution from interactions to E h (and Yj K ) at this 
level, and the corresponding Green functions are as- 
sessed by substituting E<j(f3 — t 4 ) into Eq. (7). To 
this end, we simplify the notation by setting 77 equal 
to unity, so that in the frequency space [15] E^~ + (u>) = 
-2iLue-^/ A n B (Lo), Ef(w) = 'y, a (lu)* = iuj e -^/ A , 
Ef (u) = 2iuje-^^ A (2n B (cj) + 1), with n B (u) the dis- 
tribution function for bosons. The upper frequency cut- 
off A is necessary only to ensure the convergence of the 
zero temperature parts of the Fourier integrals Hd(t) — 
J E d (a ) )(^/2 7 r)e-^ t . 

The equation of motion for the corresponding retarded 
Green function 

J-G£(ti, t 2 ) + e(p - e*Eti)G£(ti, t 2 ) = 5(h - t 2 ) (9) 
has a simple solution 

G$(h, t 2 ) = 9(h - t 2 ) cxp j- jf 2 e(p - e*Er)dr| . (10) 

In G A , the order of the time arguments is reversed in 
the step function 9(x). Consequently, after Eqs. (7), and 
(10) are combined, our problem is, in principle, solved. 



2 



To make contact with the case E — 0, we note that in 
the absence of a field, the frequency Fourier transform of 
Eq. (10) reduces to the expected result (e(p) — i-qu;)^ 1 . 

When the field is non-zero, two distinct regimes 
emerge. To obtain an estimate of the crossover field, 
we determine at which value of the field can the expo- 
nential in G R be expanded. We rewrite e(p) in the ex- 
ponential of Eq. (10) as p 2 + 2p x e*Er + m 2 . Upon in- 
tegrating the resultant exponential in Eq. (10), we find 
that the electric-field dependence can be removed from 
the exponential provided that p x e*Et 2 <C 1. Since both 
momenta and time are integrated over, this criterion is 
never satisfied. Hence, strictly speaking, the electric field 
effects are always non-perturbative. Nonetheless, if we 
focus on the characteristic timescale, t oc 1/m 2 and mo- 
mentum scale, p x oc m, we find that e*E/m 3 << 1 or 
equivalently, (e* E) 2 / 3 /m 2 = Q <C 1 determines the small 
and large field regimes. It is this criterion, rather than 
1/T 1 /* > t Ei which will serve to determine when non- 
linear transport according to Eq.(l) holds. 

To obtain the conductivity, we need first calculate 
G '". We perform the integration in the exponent of 
Eq. (10) and switch to the Wigner coordinates [14,13]: 
t = t\—t<2. and u = (ti + t 2 )/2. The Green function from 
which the current is obtained reduces to 



/>oo l>ZU 

-iG p + (t,t) = du dtM(p,e*E;,u,t)K(t) (11 

JO J-2u 

where K(t) = and 

M(p, e*E; u, t) = cxp { - (2e(p)u - 2 Px e*E{2u 2 + t 2 /2) 
(e*E) 2 



+- 



-(u 6 + 3ut 2 /2) 



(12) 



and the kernel K{t) is obtained as a result of the appro- 
priate integration over oj with the cutoff A necessary to 
regularize the T — part. The result is 



1 A~ 2 - t 2 



7T (A- 2 +i 2 ) 2 
= K (t) + K T (t) 



1 

+ - 

TT 



1 

I 2 



tt 2 T 2 



sinli (wTt) 



(13) 



in which the second temperature dependent part Kx{t) ~ 
ttT 2 /3 as Tt < 1. Similarly, because G K = G-+ + G+-, 
we can also express the inverse correlation length 



m 



p r-oo t>2u 

= S + U du dtM(p,e*E;u,t)K(t) (14) 

JpJo J-2u 



in terms of the same functions appearing in G h . 

Our problem has been reduced then to the compu- 
tation of two types of coupled integrals in Eqs. (14) 
and (11). In the T — > limit, the main difficulty 
is to obtain the cutoff-free expressions for m and J, 
since the physical results should obey the universal scal- 
ing and correspond to A -1 = 0. In this regard, we 



have found the following procedure most helpful: 1) 
Regularize the zero-temperature part (that is, the inte- 
grals over the product M(p, e*E; u, t)Ko(t)) by adding 
and subtracting M(p, e*E; u, 0) and M(p, 0;u, 0) to 
M(p,e*E;u 7 t)Ko(t). In the current, the integral over 
p in the part involving M(p, 0; u, 0) vanishes identically. 
In Eq. (11) similar contribution gives rise to the fre- 
quency and momentum dependent renormalizations of 
S, as well as the term m 2 ln(l/m 2 ), ensuring the valid- 
ity of the subsequently used logarithmic accuracy. 2) In 
the part containing M(p, e*E; u, 0) — M(p, 0;u, 0), per- 
form straightforwardly the integration over t in K (t) to 
yield 2/(wr). 3) In the part containing M(p, e*E; u, t) — 
M(p,e*E;u,0), set A -1 = as this term is com- 
pletely convergent. 4) In the terms containing the afore- 
mentioned differences, we introduce the change of vari- 
ables p x — > p x - (e*E)u/2(l + t 2 /Au 2 ), and perform the 
Gaussian integration over momenta Px and Py . Similar 
integration is performed in the part containing Kt- 5) 
Finally, change variables to y = t/2u and z — u(e*E) 2 ^ 3 . 
These transformations lead to the appearance of the aux- 
illiary functions 

2fn 2 z 

M 2 (z) = exp{- (e ^ )2/3 }, 

Mx(z, y) = M 2 (z) exp{-z 3 (l/6 + y 2 - y 4 /2)}, (15) 

and K T (z,y) = K T (t = 2z(e*E) 4 ^ 3 y). The expressions 
determining the current reduces now to quadrature: 



+ 



47T 



(e*E) i / 3 



I zdz 
Jo Jo 



dyk T (z,y)f(y)M 1 (z,y)\ (16) 



with f(y) = (1 + y 2 ). Similarly, within logarithmic accu- 
racy, we have 



m 2 ln(l/m 2 ) = A + (e*£) 2 / 3 j£° J ^ [M 

r°° dz 

-Mi{z,y)\- j o ^[M 2 (z)-M 1 (z,0)] 



2tt 



(e*£) 4 /3 



/ dz 
Jo Jo 



dyK T {z,y)Mi(z,y)\ , 



(z,0) 



(17) 



where A is the renormalized distance to the quantum 
critical point. 

Consider now the quantum-disordered regime, in 
which T is the smallest parameter. At low electric fields, 
Q = (e*E) 2 / 3 /A <C 1, we expand Eqs. (17) and (16) in 
powers of Q to obtain 



m 2 = 



and 



A (e*E) 2 tt 2 T 2 



+ 



+ 



ln(l/A) 12(A 2 /ln(l/A)) 3 A/ln(l/A) 



(18) 
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J = 



4e 



2 r 



ttT 2 
9m 4 



(e*£) 2 )l 



157rm e 



E = a(E)E, Q«l. (19) 



for the inverse correlation length and the non-linear 
current-voltage response. The first term matches iden- 
tically with the value of the linear current response cal- 
culated earlier [16] for the quantum disordered regime. 
Most importantly, the current scales as the third power 
of the voltage in the quantum disordered regime. This 
result clearly could not have been obtained within linear 
response theory. This non-linear response gives rise to a 
non-zero conductivity (albeit non-linear) even at T = 
on the putative insulating side of the transition. Such 
a term arises entirely from the zero-temperature part of 
the integrals in Eq. (16). As such a term will always be 
present regardless of which universality class describes 
the quantum phase transition, we predict that the T = 
non-linear conductivity will always remain non-zero on 
the disordered side of the transition. This conclusion is 
borne out experimentally by the series of measurements 
of dl/dV in insulator-superconductor [3,4], quantum Hall 
to insulator [5], and insulator to metal transitions [6]. 
Further, the form of the non-linearity, E 3 , is non-trivial 
and is also in agreement with the pronounced curvature 
in the I — V curves on the insulating side of the dis- 
sipative insulator-superconductor transition [4]. To ex- 
plore the large field regime from Eq. (16), we note that 
when Q 1, e*E is the only parameter that deter- 
mines to. Anticipating the logarithmic smallness, we set 
M 2 (z) w 1, and after performing first the integral over z, 
we obtain 



to = 



where 
T = 



62/3 /2 
12 I 3 



T(e*E) 2 / 3 
ln(l/(e*.E) 2 / 3 )' 

r(i + 62/ 2 -3y 4 ) 1/3 -/(2/) 



y 



= 0.116525. 



(20) 



dy 
(21) 



Here, as before, f(y) = 1 + y 2 . In the same limit, the 
current reduces to a single integral of the form, 

2e* e*E f 1 dy 2 4 

^i^^X ^mHl + ey Sy) 

e 2 

= 0.46— E Q > 1 (22) 

which proves that dl/dV approaches a universal con- 
stant, which of course depends on the universality class, 
as E — > oo. At present, the most we can conclude is that 
the constant is of order e 2 /h. This is expected to be of 
order e 2 /h as it should be on the order of the normal 
sheet resistance. Hence, we have developed a formalism 
which is capable of describing the experimental crossover 
from the non-linear to the linear regime at high field. 



The crossover to the linear regime has a fundamental 
origin. Because at mean- field v = 1/2 and z = 2, we can 
rewrite Eq. (20) as to 2 oc 1/ {i 1 ^ \n{^ v )). Hence, at 
large fields, the correlation length is cutoff by the elec- 
tric field length scale, to logarithmic accuracy. That is, 
the only length scale in the large field limit is Ie- As it 
is the product (uiIe) v that enters the non-linear conduc- 
tivity, the electric field dependence naturally drops out 
and the conductivity approaches a universal constant. In 
the opposite regime, the correlation length and the elec- 
tric field length scale are distinct as Eq. (18) indicates. 
The non-linearity arises entirely from the quantum fluc- 
tuations on the length scale £. In the non-linear regime, 
it is tempting to rewrite the current, Eq. (19), in terms 
of an effective temperature, T 2 ff = T 2 + 0.06(e*£ym) 2 . 
The current then simplifies to 



47re* 
9h 



tLe. 



(23) 



For electrons localized in the band tails of semiconduc- 
tors, the concept of an electric field-dependent effective 
temperature has been used extensively [18,19]. Electrons 
moving against the electric field are accelerated on a 
length scale set by the localization length. In the quan- 
tum disordered regime, it is the correlation length that 
plays the role of the localization length. Consequently, 
it is not surprising that the bosonic excitations can be 
described by an effective E-field dependent effective tem- 
perature. 
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